The Theory of p-n Junctions in Semiconductors and p-n 
Junction Transistors 

By W. SHOCKLEY 

In a single crystal of semiconductor the impurity concentration may vary 
from //-type to »-type producing a mechanically continuous rectifying junction. 
The theory of potential distribution and rectification for p-n junctions is developed 
with emphasis on germanium. The currents across the junction are carried by 
the diffusion of holes in n-type material and electrons in />-type material, re- 
sulting in an admittance for a simple case varying as (1 + mot,,)" 2 where r p is the 
lifetime of a hole in the H-region. Contact potentials across p-n junctions, carry- 
ing no current, may develop when hole or electron injection occurs. The principles 
and theory of a p-n-p transistor are described. 
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1. Introduction 

AS IS well known, silicon and germanium may be either »-type or 
//-type semiconductors, depending on which of the concentrations 
Nd of donors or A 7 ",, of acceptors, is the larger. If, in a single sample, there 
is a transition from one type to the other, a rectifying photosensitive p-n 
junction is formed. 1 The theory of such junctions is in contrast to those 

1 For a review of work on silicon and germanium during the war see H. C. Torrey and 
C. A. Whitmer, Crystal Rectifiers, McGraw-Hill Book Company, Inc., New York (1948). 
P-n junctions were investigated before the war at Bell Telephone Laboratories by R. S. 
Ohl. Work on p-n junctions in germanium has been published by the group at Purdue 
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of ordinary rectifying junctions because, on both sides of the junction, 
both electron flow and hole flow must be considered. In fact, a major 
portion of the hole current may persist into the w-type region and vice- 
versa. In later sections we show how this feature has a number of inter- 
esting consequences, which we shall describe briefly in this introduction. 

A p-n junction may act as an emitter in the transistor sense, since it can 
inject hole current into w-type material. The a-c. impedance of a p-n junc- 
tion may exhibit a frequency dependence characterized by this diffusion 
of holes and of electrons. For high frequencies the admittance varies ap- 
proximately as (Vw) 1 ' 2 and has comparable real and imaginary parts. When 
a p-n junction makes contact to a piece of n-type material containing a high 
concentration of injected holes, it acts like a semipermeable membrane and 
tends to come to a potential which corresponds to the hole concentration. 

Although some results can be derived which are valid for all p-n junctions, 
the diversity of possible situations is so great and the solution of the equa- 
tions so involved that it is necessary to illustrate them by using a number 
of special cases as examples. In general we shall consider cases in which the 
semiconductor may be classified into three parts, as shown in Fig. 1. The 
meaning of the transition region will become clearer in later sections; in 
general it extends far enough to either side of the point at which N d — N a = 
so that the value of | N d - N a \ at its boundaries is not much smaller than 
in the low resistance parts of the specimen. As stated above, appreciable 
hole currents may flow into the w-region beyond the transition region. For 
this reason, the rectification process is not restricted to the transition region 
alone. We shall use the word junction to include all the material near the 
transition region in which significant contributions to the rectification 
process occur. It has been found that various techniques may be employed 
to make nonrectifying metallic contacts to the germanium; when this is 
properly done, the resistance measured between the metal terminals in a 
suitably proportioned specimen is due almost entirely to the rectifying 
junction up to current densities of 10 -1 amp/cm 2 . 

directed by K. Lark-Horovitz: S. Benzer, Phvs. Rev. 72, 1267 (1947); M. Becker and 
H. Y. Fan, Pliys. Rev. 75, 1631 (1949); and H. Y. Fan, Phys. Rev. 75. 1631 (1949). Similar 
junctions occur in lead sulfide according to L. Sosnowski, J. Starkiewicz and O. Simpson, 
Nature 159, 818 (1947), L. Sosnowski, Phys. Rev. 72, 641 (1947), and L. Sosnowski, B. 
W Socle and J. Starkiewicz, Nature 160. 471 (1947). The theory described here has been 
discussed in connection with photoelectric effects in p-n junctions by F. S. Grmcher.. 
Meeting of the American Physical Socielv, Cleveland, March 10-12, 1949 and by W. 
Shocklev, G. L. Pearson and M. Sparks. Phys. Rev. 76, 180 (1949). For a general review 
of conductivity in p- and w-tvr>e siHcon see G. L. Pearson and J. Bardeen, Phys. Rev. 75, 
865 (1949) and |. H. Scaff, H. C. Theuerer and E. E. Schumacher, Jl. of Metals, 185, 
383 (1949)'and YV. G. Pfann and J. H. Scaff, Jl. of Metals, 185, 389 (1949). The latter 
two papers also discuss photo-voltaic barriers. The most recent and thorough theory for 
frequency effects in metal semiconductor rectifiers is given elsewhere in this issue Q. 
Bardeen, Bell Sys. Tech. Jl., July 1949). 
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Even for distributions of impurities as simple as those shown in part (b) 
there are two distinctly different types of behavior of the electrostatic po- 
tential in the transition region, each of which may be either rectifying or 
nonrectifying. The requirement that the junction be rectifying can be stated 
in terms of the current distribution, certain cases of which are shown in (c). 
The total current, from left to right, is /, the hole and electron currents being 




DISTANCE THROUGH SAMPLE , X >■ 

Fig. 1 — -The p-n junction. 

(a) Schematic view of specimen, showing non-rectifying end contacts and convention 
for polarities of current and voltage. 

(b) Distribution of donors and acceptors. 

(c) Three possible current distributions. 

I p and /„ , with / = I p -\- I„. Well away from the junction in the p-type 
material, substantially all of the current is carried by holes and I p = /; 
similarly, deep in the «-type material I„ = I and I p = 0. In general in a 
nonrectifying junction, the hole current does not penetrate the w-type ma- 
terial appreciably whereas in the rectifying junction it does. Under some 
conditions the major flow across the junction will consist of holes; such 
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cases are advantageous as emitters in transistor applications using «-type 
material for the base. 

Where the hole current flows in relatively low resistance n-type material, 
it is governed by the diffusion equation and the concentration falls off as 
exp(-x/L p ) where L p is the diffusion length: 

L= VdT p . 

Here D is the diffusion constant for *holes and r p their mean lifetime. The 
lifetime may be controlled either by surface recombination 2 or volume re- 
combination. Surface recombination is important if the specimen "has a 
narrow cross-section. 

Under a-c. conditions, the diffusion current acquires a reactive component 
corresponding to a capacity. In addition, a capacitative current is required 
to produce the changing potential distribution in the transition region 

itself. 

In the following sections we shall consider the behavior of the junction 
analytically, treating first the potential distribution in the transition region 
and the charges required change the voltage across it in a pseudo-equilibrium 
case. We shall then consider d-c. rectification and a-c. admittance. 

2. Potential Distribution and Capacity of Transition Region 

2.1 Introduction and Definitions 

We shall suppose in this treatment that all donors and acceptors are 
ionized (a good approximation for Ge at room temperature) so that we have 
to deal with four densities as follows: 

n = density of electrons in conduction band 
p = density of holes in valence-bond band 
Na = density of donors 
N a = density of acceptors 
The total charge density is 

p'- q (p - n + N d - N tt ), (2.1) 

where q is the electronic charge. We shall measure electrostatic potential 
^ in the crystal, as shown in Fig. 2, from such a point, approximately mid- 
way in the energy gap, that if the Fermi level <p is equal to \p, the concentra- 
tions of holes and electrons are equal to the concentration »,- = pi char- 

' H. Suhl and W. Shockley Phys. Rev. 75 1617 (1949). 

3 A difference in effective masses for holes and electrons will cause a shift of \f> from the 
midpoint between the bands, 
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Fig. 2 — Electrostatic potential ip, Fermi level <p and quasi Fermi levels <p p and y>„. 
(In order to show electrostatic potential and energies on the same ordinates, the ener- 
gies of holes, which are minus the energies of electrons, are plotted upwards in the figures 
in this paper.) 

acteristic of a pure sample. For an impurity semi-conductor we shall have, 
as shown in (b), 



where q is the electronic charge. Accordingly, 

P = q\N d - N + Irii sinh [qfa - +)/kT\}. 



(a) 
(b) 



(2.2) 



(2.3) 



When the hole and electron concentrations do not have their equilibrium 
values, because of hole or electron injection or production of hole-electron 
pairs by light, etc., it is advantageous to define two non-equilibrium quasi 
Fermi levels <p v and tp n by the equations 

(a) , X 

(2.4) 
(b) 

as indicated in Fig. 2 (c). In terms of tp p and <p„ , the hole and electron cur- 
rents take the simple forms: 

I„ = -q[DVp + ppty] = -qnpV<p p (2.5) 

/„ = bq[DVn - /i»V^] - -qbnnVip n (2.6) 



p = me" 1 *-*""* 



440 BELL SYSTEM TECHNICAL JOURNAL 

where the mobility p and diffusion constant D for holes are related by Ein- 
stein's equation 

H = qD/kT (2.7) 

and b is the ratio of electron mobility to hole mobility. 4 

Under equilibrium conditions <p p = <p n = <P where <p is independent of 
position. Under those conditions, I p and /„ are both zero according to equa- 
tions (2.5) and (2.6). The electrostatic potential ^, however, will not in 
general be constant and there will be unbalanced charge densities throughout 
the semiconductor. We shall consider the nature of the conditions which 
determine ty for a general case and will later treat in detail the behavior of 
\f/ for p-n junctions. 

For equilibrium conditions, there is no loss in generality in setting <p 
arbitrarily equal to zero. The charge density expression (2.3) may then be 
rewritten as 

P = pa — Pi sinh u (2.8) 

where 

U = qt/kT, pi = 2mq, p d = q(N d - N a ) (2.9) 

In equation (2.8) p d and u and, consequently, p may be functions of position. 
The potential \f/ must satisfy Poisson's equation which leads to the equation 

VV = -W< (2.10) 

where k is the dielectric constant, (2.10) can be rewritten as 

V > u J^Uhu- P -A (2.11) 

fc£7\ \ pi/ 

What this equation requires in physical terms is that the electrostatic po- 
tential produces through (2.8) just such a total charge density p that this 
charge density, when used in Poisson's Equation (2.10), in turn produces 
ip. It seems intuitively evident that the equation for u will always have a 
physically meaningful solution; no matter how the charge density p d due 
to the impurities varies with position, the holes and electrons should be 
able to distribute themselves so that equilibrium is produced. For a 
one-dimensional case, it is not difficult to prove that a unique solution exists 
for u{x) for any p d (x) (Appendix VII). 

4 We prefer b in comparison to c for this ratio since c for the speed of light also occurs in 
formulae involving b. 
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The coefficient in (2.11) has the dimensions of (length) -2 leading us to 
define a quantity 

L D = VKkT/Airqpi = VKkT/Sirq 2 /^ 

= 2.1 X 10~ 3 cm for Si with k = 12.5, 8 «, = 2 X 10 10 cm -3 (2.12) 

= 6.8 X 10~ 5 cm for Ge with k = 19, 6 m = 3 X 10 13 cm -3 

where the subscript D for Debye emphasizes the similarity of L D to the char- 
acteristic length in the Debye-Hiickel theory of strong electrolytes. The 
meaning of the Debye length is apparent from the behavior of the solution 
in a region where pa is constant, and u differs only slightly from the value 
wo which gives neutrality, with p,- sinh w = Pd • Under these conditions, 

— = (LI? cosh Mo)(m - tt ) (2.13) 

so that u — i<o varies as exp (± .r\/cosh tio/L D ). In general, we shall be in- 
terested in cases in which the deviation of u from «o decays to a small value 
in one direction. It is evident that the distance required to reduce the devia- 
tion to l/e is L D / Vcoshl/o . If only small variations in p d occur within a dis- 
tance L D I \/cosh Wo , then the semiconductor will be substantially neutral. 
However, if a large variation of pa occurs in this distance, a region of local 
space charge will occur. These two cases are illustrated in connection with 
the potential distribution in a p-n junction. 

2.2 Potential Distribution in the Transition Region 1 

We shall discuss the case shown in Fig. 1 for which the charge density 
due to donors and acceptors is given by 

N d - N a = ax (2.14) 

This relationship defines a characteristic length L a given by 

L a = ru/a (2.15) 

If L a ~2> L D , the condition of electrical neutrality is fulfilled (Appendix VII) 
and u satisfies the equation 

sinh u = pa/pi — ax/2>ii = x/2L a 

6 J. F. Mullaney, Phys. Rev. 66, 326 (1944). 

6 H. B. Briggs and W. H. Brattain, Pliys. Rev., 75, 1705 (1949). 

7 Potential distributions in rectifying junctions between semiconductors and metals 
have been discussed by many authors, in particular N. F. Mott, Proc. Roy. Soc. 171A, 
27 (1939) and VV. Schottky Zeits.f. Physik 113, 367 (1939) 118, 539 (1942) and elsewhere. 
A summary in English of Schottky's papers is given by J. Joffe, Electrical Communications 
22, 217 (1945). All such theories are in principle similar in involving the solution of equa- 
tions like (2.11). See, for example, H. Y. Fan, Phys. Rev. 62, 388 (1942). 
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On the other hand, if L D » Z« , a large change in impurity concentration 
occurs near x = without compensating electron and hole densities oc- 
curring. Mathematically, we find that (2.11) can be expressed in the form 

— =-i(-y + sinh«) (2.16) 

dy- A 2 

and 

y = x/2L a , K = U/lLa (2.17) 

In Appendix VII, it is verified that the appropriate solution for K « 1 is that 
giving local neutrality, u = sinh -1 y; while for K » 1, there is space charge 
as described below. 

For L D » L a , or A' » 1, there is a space charge layer in which iXj — N a 
is uncompensated. To a first approximation, we can neglect the electron and 
hole space charge in the layer and obtain, by integrating twice, 

_ _2w^ + (218) 

OK 

where we have chosen the zero of potential as the value at x = 0, a condi- 
tion required by the symmetry between +.r and —x of (2.14). Although 
the potential rise is steep in the layer, d\J//dx should be small at the point 
Xm where the neutral «-type material begins. As an approximation we set 
dift/dx = at .v = x m : 

d ±= _ 2 W" + fl2 = o ; (2.19) 

dx K 

this leads to a value for o 2 which may be inserted in (2.18) to evaluate yf/ 

at x m : 

4Trqax m 4-rrq , x z 4irq 3 , n „„■> 

where n m = ax m is the density of electrons required to neutralize Nd — N a = 
ax m at the edge of the space-charge layer. This value of n,„ must corre- 
spond to that associated with \p m by (2.2) 

n m = tii e 

We thus have two equations relating rp m and //„, and the parameter "a." 
To solve them we plot ln^ m versus In n„, as shown in Fig. 3. On this figure the 
relationship 



_ ., . „9*m/*T (2 21) 



47ro n m 



3 

3k a? 



3 

= 3.18 X 10~ 8 ^ volts for Ge 
a 2 



= 4.83 X 1(T 8 ^ volts for Si (2.22) 
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becomes a family of straight lines with "a" as a parameter. (Only a = 10 
cm -4 is shown for Si, all the other lines being for Ge.) The half thickness 
x m (= ii m /a) of the space-charge region is also shown. Solutions are obtained 
when these lines cross the curves n m = »< exp (q\f/ m /kT), which are shown 
for room temperature. The condition that the intersection lie well to the 
right on the curve is equivalent to K » 1. For two Si samples cut from a 
melt, a was determined from measurements of conductivity and was 
about 10 15 to 10 16 cm -4 . For these, the space change region has a half-width 
x m of more than 10~ 4 cm. For other temperatures, the curves can be ap- 
propriately translated. 9 
In Fig. 4(a) we show the limiting potential shapes: 

ax = 2m sinh % for K « 1 (2.23) 

* = 6W2)(- (*/*«)' + 3(x/x m )) for K » 1 (2.24) 

In Fig. 4(b) the charge densities are shown. For the space-charge case, 
| Nd — A r | is greater than n or p. For a higher potential rise, i.e. larger 
}f/ m , the discrepancy would be greater and Na — N a would be unneutralized 
except near x m . 

2.3 The Transition-Region Capacity 

When the voltage across the junction is changing, a flow of holes and 
electrons is required to alter the space charge in the transition region. We 
shall calculate the charge distribution in the transition region with the aid 
of a pseudo-equilibrium model in which the following processes are imagined 
to be prevented: (1) hole and electron recombination, (2) electron flow across 
the />-region contact at .r a (Fig. 1), (3) hole flow across the «-region boundary 
at x b . Under these conditions holes which flow in across .r a must remain in 
the specimen. If a potential 8<p is applied at the p end, then holes will flow 
into the specimen until <p p has increased by 8<p so that the holes insiie are 
in equilibrium with the contact which applies the potential. Since the speci- 
men as a whole remains neutral, an equal electron flow will occur at x b . 
When the specimen arrives at its pseudo-equilibrium steady-state, the 
potential distribution will be modified in the transition region and the num- 
ber of holes in this region will be different from the number present under 
conditions of true thermal equilibrium. The added number of holes is pro- 
portional to 5<? for small values of b\p and thus acts like the charge on a con- 
denser. Our problem in this section is to calculate how this charge depends 

8 Unpublished data of W. H. Brat tain and G. L. Pearson. 

9 The effect of unionized donors and acceptors can also be included by letting tu include 
the properly weighted donor states and />,-, the acceptor states. 
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upon 8<p for various types of transition regions and to express the result as 
a capacity. 

The justification for this pseudo-equilibrium treatment is as follows: 
Under actual a-c. conditions the potential drop in the p- and n-regions them- 
selves are small because of their high conductivity so that most of the po- 
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Fig. 4 — Electrostatic potential and densities for p-n junctions. 

tential drop occurs across the transition region. On the p-side of the transi- 
tion region a large supply of holes is available to modify the potential and 
the fact that a current is flowing across the junction disturbs their concen- 
tration negligibly; the electrons on the «-side are similarly situated. Hence 
the distribution of holes and electrons in the transition region will be much 
the same as for the pseudo-equilibrium case. The question of how the hole 
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current required to change the potential distribution in the transition region 
is related to other hole currents is discussed in Section 4.1. 

Under our assumptions, after the voltage 8ip is applied, a steady state is 
reached involving no current hence V<p p = V(p n = 0. Consequently, both 
(p p and (p n are constant and 

<P P — <p n = 8<P (2.25) 

since the holes are being supplied from a source at a potential 5^ higher than 
for the electrons. 
We shall then have 

p = me* 19 *-*""* = nie* (n '- mT (2.26) 

n = n i e q{ *- pn),hT = n 1 e«*-* i) ' hT (2.27) 
where 

<Pi = (<P P + <Pn)/2, <fp = tpi + Sip/2, ip n = tpi — titp/2 (2.28) 

and 

(2.29) 



fl 5v>/2fcr 
iii = Hie 



Thus the effect of applying the potential 8<p in the pseudo-equilibrium case 
is equivalent to changing «,• to Mi just'as if the energy gap had been reduced 

by qbip. 

In the />-region, n « p and so that p = — ax is a good approximation. 
Similarly, in the w-region, we set n = ax. Hence we have in the /(-region 

^ = n + fy/2) - (kT/g) In (-ax/m) (2.30) 

and in the ^-region 

* - ifii - (fr/2) + (kT/q) In (ax/m). (2.31) 

Hence the effect of dip is to shift \j/ in the ^-region upwards by 8<p compared 
to ^ in the w-region. This is an example of the general result that \p — <p p 
tends to remain constant at a given point in the />-region no matter what dis- 
turbances occur and \j/ - <p n tends to remain constant in the ra-region. 

The Capacity for the Neutral Case K « 1 

For the neutral case, we calculate the total number of holes, P, between 
x c and x b as a function of &p. The charge of these holes is qP and the effective 
capacity is q dP/d dip. As explained above, we are really interested in the 
change in number of holes in the transition region. However, the value of P 
is relatively insensitive to the location of the limits x a and x b so long as they 
lie in regions where the conductivity approaches the maximum values in the 
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p- and ^-regions. In the following calculations, we shall consider a unit area 
of the junction so that values of P and of capacity are on a unit area bases. 
The value of P is obtained by integrating p dx making use of the neu- 
trality condition to establish the functional relationship between p and x. 
The neutrality condition can be written as 

ax = 2wi sinh Z, = 2m sinh u (2.32) 

where // = qty - ipi)/kT and 

p-**tor*'* r mn l n (2.33) 

n = >he +u (2.34) 

so that the value of P can be obtained by changing variables from x to u: 

p*b put, 

P = / p dx = I p(2ni/a) cosh u du 

r . b (2.35) 

= W/fl) J u [1 + e~ 2u ] du = (»!/«) [«6 - M. + (r*- - e- 2 "»)/2] 

For the cases of practical interest, the value of p at x = x a , denoted by 
p a , and the value of n at * = x b , denoted by n b , will both be large compared 
to Hi . Consequently, we conclude that 

u a = —\n(p a /m) and u b = In n b /n x 

are both larger than unity in absolute value but probably less than twenty 
for a reasonable variation of impurity between x a and x b . (For example for 
a change in potential of 0.2 volts such as would occur between p- and w-type 
germanium, u a and u b would each be about 4 in magnitude.) Hence we ob- 
tain for P, 

P = (nl/2a)(2(ti b - u a ) + {p a / m f - (;/ x /« 6 ) 2 ) 

££ pi/ 2a + (th/o)(ut, - u a ) (2.36) 

where we have neglected {n x /n b f which is «1 and the negligible compared 
to u b — u a . The term pa/ 2a is simply the integrated acceptor-minus-donor 
density in the p-reg\on, as may be seen as follows: 

j (A r a - iVd) dx = jT (-ax) dx = axl/2 = p\/2a. (2.37) 

The second term in (2.36) is essentially the sum of the holes of the right 
of x = plus the electrons to the left of x = 0, whose charge is also com- 
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psnsated by holes. The total number of holes can be expressed in terms of 
5<p through the dependence of »i on 8<p. The second term is thus 

(n]/a)[\n(n b /th) + In (p a /ni)] 

= {n\/a)e 9i9,kT • [\n(n b p a /n)) - q8?/kT] (2.33) 

Hence for a small change d8tp in 8<p, the change in charge dQ= qdP and the 
capacity C are given by 



2 2 



m K\qWb- ft, ~ W/*± - Ug (241) 



C = |2 - » 2! [In^^iS - (qS-f/kT) - 1]. (2.39) 

This capacity can be reexpressed in terms of the difference in \p between 
x a and x b : When 8<p = 0, corresponding to the thermal equilibrium case, 
we have 

p a n b = nW l * b -* a) ' kT (2-40) 

Using this together with the definitions of L D and L a we obtain 

ifoCfc - + a - 8f )/kT - 1] e** ,hv 

In this expression ^ and fa are the potentials when 8<p = 0; so that 

^6 — (<Pa + M 

is thus the increase in potential in going from Xa to x b when Sip is applied. 

For thermal equilibrium, 8>p = and, as discussed above, the term in 

fo _ fa w ju be about 10. Hence, using the definition K = L D /2L a , we have 

C ££ k/4tt(4^L d /10) (2.42) 

For if « 1, the case for which this formula is valid, C will be the capacity 
of a condenser whose dielectric layer is much less than L D thick. 

Capacity for Space Charge Case, K^> 1 

As discussed in connection with (2.30) and (2.31), the applisd potential 
8<p reduces the increase (= 2\f/ m ) in rp between the ^-region and the »-region 
by &p/2 on each side of x = 0. This is accomplished by a narrowing of the 
space charge layer by 8x m on each side where (according to (2.20)) 

ty m = -V/2 = brqaxUxm/it (2.43) 

The decrease in width 8x m brings with it an increase in number of holes 
-ax Sx m per unit area of the junction on the p-side and an equal number of 
electrons on the »-side. Thus a charge of holes psr unit area of 52 = — qax m 8x m 
must flow in from the left. The capacity per unit area is, therefore, 

C = 8Q/8<p = qax m 8x m /8<p - k/^t2x„ (2.44) 
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corresponding to a condenser of thickness 2x m . It is evident that formula 
(2.44) will hold for a small change d8ip superimposed on a large biis 8<p pro- 
vided that 2x m is the thickness of the space charge region under the condi- 
tions when 8<p is applied. If \f/„ is the value of \f/ for 8<p = 0, then ^ m = 
V'n.o — Sip/2-, and C will vary as 

C = K[4*qa/3Kty„. - V/2)] 1/3 /8tt (2.45) 

so that 1/C 3 should plot as a straight line versus Sip with slope 

- (8irA) 8 (3«/8irsa) = - i-^ . (2.46) 

a'qa 

In addition to the holes which flow to account for the change in \f/ m , the 
concentration of holes in the «-region will be increased by a factor 
exp((j8tp/kT). However, this increase does not lie in the transition region; 
we shall consider it later, in Section 4, in connection with a-c. admittance. 

Comparison of the Two Capacities 

It is instructive to compare the two capacities just derived. We suppose 
that for one value of w, we have K » 1 so that the space charge solution is 
good. For this case we choose x a = —x m and xi, = -\-x m so as to bound the 
space charge layer. We then imagine w, to be increased, either by raising 
the temperature or by applying a potential difference 8ip. The capacity then 
changes from 

C B p. c hg. = fc/87T.V m tO Cneut. = 5k/8tKL D (2.47) 

(i.e., from (2.44) to (2.42)) so that the ratio is 

C neut. oX m 

c~^r = kT (2 - 48) 

For K < 1, this ratio is large, both because of K in the denominator and 
because x m > L a so that x m /L D > L a /L D = 1/2 K. 

In Section 4.4 we shall compare these capacities with that due to dLIusion 
of holes and electrons beyond the transition region. 

2.4 The Abrupt Transition 

For completeness we shall consider the case in which the impurity con- 
centration changes abruptly from p p to n„ at x = 0. For this case the po- 
tential in the space-charge layer will be of the parabolic type discussed by 
Schottky, the potentials varying as 

\p = (2ir/K)q p p (x - x p f + constant, x < (2.49) 

f = — {2r/K)q n„(x — x n ) 2 + constant, x > (2.50) 
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where *, < and x n > are the ends of the space-charge layer in the p- 
and w-regions. The gradient of potential at x = must be equal for the two 
layers leading to 

—ppXp = n n x n (2.51) 

so that if the total width of the space charge layers is W = x n - x p , it 
follows that 

*, = -n n W/(n n + p p ) and .r n = p p W/(n n + p p ). (2.52) 

The potential difference across the layer, which is \p b — fa is 
fa - fa = (2t ? A)(/>p4 + »»*») = [2irq p P n,,/K(p p + n n )]W 2 , (2.53) 
If p p » n„ this reduces to 

fa-fa = lirq ,i„W 2 /k (2.54) 

the formula given by Schottky, which should be appreciable in this case, 
for which all the voltage drop occurs in the w-region. 
The capacity for the abrupt transition will be 

C = k/AttW (2.55) 

where W is obtained by solving (2.53). For this case (1/C) 2 should plot as a 
straight line versus \p b — ^. : 

— = \8w(p p + n„)/Kqp p nMb ~ to). (2.56) 

3. General Conclusions Concerning the Junction Characteristic 

In this section we shall consider direct current flow through the junction 
and shall derive the results quoted in Fig. 1 relating the current distribution 
to the characteristics of the junction. We shall suppose that holes and elec- 
trons are thermally generated in pairs at a rate g and recombine at a rate 
rnp so that the net rate of generation per unit volume is 

(net rate of generation) = g — rnp, 

which vanishes at equilibrum. Obviously, g = rn* . If relatively small con- 
centrations 8p and 5n of holes and electrons are present in excess of the 
equilibrium values, the net rate of generation is 

S'p = 8h = g - r(n + Sii)(p + Sp) = -rnSp - rpSn (3.1) 

This is equivalent to saying that excess holes in an «-type semiconductor, 
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and excess electrons in a p-type semiconductor, respectively, have lifetimes 
t p and t„ given by 

8'p = —8p/T p = —rn8p or t p = l/rn = p/g (3.2) 

and 

8h = —8h/t„ = —rp8n or t„ = l/rp = n/g. (3.3) 

We shall have occasion to use this interpretation later. (We later consider 
the modifications required when surface recombination occurs, Section 4.2, 
Appendix V, and the effect of a localized region of high recombination rate, 
Section 4.6, Appendix III.) 

In principle, the steady-state solution can be obtained in terms of the 
three potentials \p, (p p and <p n . These must satisfy three simultaneous or- 
dinary differential equations, which we shall derive. As discussed in Section 
2, we consider all donors and acceptors to be ionized so that Poisson's equa- 
tion becomes 

*+ = J^l ( ax + me^-*" kT - me**-*^*) (3.4) 

dx 2 k 

an equation in which the unknowns are the three functions <p p , <p„ and yp. 
The total current density, from left to right, is 



/ = /, + /„ = -qn 



[>% + *%]■ ^ 



The elimination of p and n by equation (2.4) results in an equation in- 
volving the three unknown functions and I. The divergence of hole current, 
equal to the net rate of generation of holes, is 



- -nut [j- ( d<fi p\ _ ^_dypdj>p d 2 <p p ~\ 
WV [_kT\dx) kTdxdx dx 2 ] 



dh 

dx "' \_k'i \dx / HI Ax dx dx* J (3 6) 

= q(g - rnp) = qg(\ - e "^-^ lkT ), 



with p in the second term given by (2.4) so that (3.6) is also an equation 
for the three unknown functions. The equation for dl„/dx can be derived 
from the last two and adds nothing new. These three equations can be used 
to solve for aip/dx , a<p p /dx and d<p„/dx in terms of lower derivatives and /. 
They thus constitute a set of equations sufficient to solve the problem pro- 
vided that physically meaningful boundary conditions are imposed. We 
shall not, however, deal directly with these equations; the main reason for 
deriving them was to show that the problem in question is, in principle, 
completely formulated. Instead of attempting to solve the equations, we 
shall discuss certain general features of the solutions for ip p and <p„ , using 
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approximate methods, and in this way bring out the essential features of 
the theory of rectification. 




DISTANCE THROUGH SAMPLE , X— *- 
Fig. 5 — Potential and current distributions for forward current in p-n junctions. 

(a) p-n junction under equilibrium conditions. 

(b) Division of current between holes and electrons. 

(c) Distribution of potentials for forward current flow showing how the potential 5? 
applied at .r a changes <p p , <p„ and <p. 

In Fig. 5 we represent a general situation which may be used to illustrate 
the nature of the resistance of the junction. Part (a) corresponds to thermal 
equilibrium and shows the potential distribution and Fermi level in ac- 
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cordance with the scheme used in Fig. 2. Part (b) shows the current dis- 
tribution for a forward current I from left to right and (c) shows 
the corresponding potential distribution and values of <p n and <p p , the total 
applied potential being 8<p. Recombination prevents the hole current from 
penetrating far into the «-region, the depth of penetration being described 
by the diffusion length L p = ~s/Dt p = s/Dp n /g, where p n is the hole con- 
centration in the w-region. The electron current similarly is limited by 
L n = \ZbDr„ = \/bDn p /g. (Diffusion lengths are evaluated for particular 
models of the junction in Section 4.) Far from the junction, therefore, the 
hole and electron concentrations have their normal values and consequently 
<P P = <p n and (f p — }f/ has its normal value. This accounts for the equal dis- 
placement 8<p for all three curves at x = x a . The curves for <p p and ip n have 
a continuous downward trend which produces the currents 

I p =-qtf*!g and I n =-qb m d ^. (3.7) 

The area between the <p p and <p„ curves has a special significance: This differ- 
ence is related to the excess rate of recombination and the integral of this 
rate over the entire specimen must be sufficient to absorb the hole current 
I P = I entering at x„ so that the entire current at Xb is carried by electrons. 
In terms of <p p — <p n and equation (3.6) we obtain 

I = I p (x a ) - I p (xb) = I -dl p 

(3.8) 



- M [ (fi 9{ **-^' kT - 1) dx. 

"In 



From (3.8) we conclude that if g is increased indefinitely for a specified 
current /, then <p p — tp n must approach zero. For this case, in which the 
rate of -recombination and generation is very high, <p p = <p n and 



and 



bifi 



I = I p + I n = -qn(p + hi) dipjdx (3.9) 



= - ! d<p p = I f dx/qy.{p + bn) = IR , (3.10) 



where Bo is simply the integral of the local resistivity corresponding to 
densities p and n. For smaller values of g, I does not divide in the ratio p'.bn 
and <p p 9* <p n and hip > IRq . 

We shall next give an approximate treatment for the case in which dpj 
(J for junction), the value of <p p — <p n at x = 0, is an appreciable fraction of 

10 A general proof that 6<p > IRo is given in Appendix I. 
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the total voltage drop. For this purpose we treat ip p — <p H as constant over a 
range of integration from x = — L n to .v = +L P obtaining 

I = gq(L n + L}\f*» m - 1] 

- UP**'** - 1| 

where 

/.= gq(L n + L p ) (3.12) 

is the current density corresponding to the total rate of generation of hole- 
electron pairs in a volume L n + L P thick. We next consider 8<p Kp + 8ip Sn 
shown in Fig. 5c, where, as the subscript R implies, these are thought of as 
resistive terms and are given by the integrals 

8<Pr p + 8<p Rn = - / <h P - I dtp n = / I v dx/qtip + / /„ dx/qubn. 

The denominators are both approximately qn(p + bn) which occurs in the 
integral for Rq . Furthermore, for most of the first range I p = I and for 
most of the second I n = I. Near x = 0, I p or /„ must be at least 1/2. Hence 
it is evident that 8cp Kp + fyitn cannot be much less than IR ■ We shall repre- 
sent it by IRi where R < 2Ri < 2R . 

In terms of Ri and /„ , the relationship between current and voltage 
becomes 



* bT / j"\ 

Sip = StpBp + S<PRn + 8<pj = Ril + — In (l + J ) ■ 



(3.13) 



This corresponds to an ideal rectifier in series with a resistance Ri . The 
junction will, therefore, be a good rectifier if the second term represents a 
much higher resistance. 

We shall compare the two resistances for the case corresponding to K « 1 . 
For this case, we have p = -ax and n = +ax except in the narrow range 
\x\ < L a = Hi/ a. The integral R can be approximated by integrating 
dx/tr for x outside of the range ±Z. a using the approximation ±a.r for 
p and n and approximating the integral from — L a to +L a by 2L a /a (in- 
trinsic). This procedure gives 



r~ Xa 2L f 

Rx = / dx/qnax + . ° > + / dx/qpbax 

(3.14) 

= — (l + Tjln(VO 
qilfli \ b/ 

where it is supposed that -x a = x b and that In (x b /L a ) is large compared to 
2/(b + 2 + l/b). The evaluation of L p and L n for use in I 8 is more involved 
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since t p and t„ are both functions of x. We shall obtain an approximate self- 
consistent diffusion length by assuming that the holes diffuse, on the average, 
to just such a depth, L p , that in uniform material of the type found at 
L p , their diffusion length would also be L p . At a depth L p , the value of 
n is aL p so that by (3.2), t p is \/raL p = n"i/gaL p . Thus we write 

L% = Dt p = DnVgaL p . (3.15) 

We can solve the equation (3.15) for L p and a similar one for L n and insert 
the results in equation (3.13). For small / this gives 

6<p/I = R t + {kT/qh) = A (l + t) In ( Xb /L a ) 

?M«. \ bj ( 3 16 ) 

+ kT/H g 2 '\DL a n,) m a + b 1 ' 3 )). 

It is seen that for g large, the second term, corresponding to the rectifying 
resistance, becomes small. For this case, as discussed above, <pp = <Pn and 
the exact integral for i?o should be used and the junction will give poor 
rectification. 

It is also instructive to consider L a as a variable. Increasing L a corre- 
sponds to making the transition from p to n more gradual. It is evident 
that varying L a changes the two terms of (3.16) in opposite directions so 
that there will be an intermediate value of L a for which the resistance of 
the junction is a minimum. As L a approaches zero, however, the second 
term should be modified: If we imagine that in the transition region the 
concentration (Nj — A r a ) varies only over a finite range, bounded by fixed 
values ;/„ and p p in the n- and />-regions, then it is clear that the limiting 
values of L p and L n should be given not by (3.15) but by -\/Dt p and y/bDr n 
where t p and t„ are evaluated in the //-region and />-region. This leads to a 
limiting value for I s , which is given in equation (4.11) of the following 
section. In the range for which (3.16) applies, however, the interesting 
result holds that widening the transition region initially decreases the re- 
sistance by furnishing a larger volume in which holes and electrons may com- 
bine or be generated. 

The condition that 8tpj dominate the resistance is that the second term 
of (3.16) be much larger than the first. This leads to the inequality 

1 « tr<SL*v ■ xr' = (fi *" ar ' (L " /L ' r (3 ' 17) 

where we have neglected various factors involving b, which are nearly unity, 
and \n(xb/L a ) (which must be about 4 for Ge since the conductivity at Xb 
is about exp(4) times the intrinsic conductivity). The quantity 

L pi = (Dm/gf 2 (3.18) 
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is the diffusion length for holes in the intrinsic region. The inequality states 
that the diffusion length must be much larger than L a . This is equivalent 
to the previous statement that the hole current must penetrate the w-region 
for the rectifier to have a good characteristic. (If a local region of high re- 
combination is present in the transition region, this result just quoted need 
not apply. See Section 4.6.) 

If the hole current penetrates deeply into the w-region and Ri is negligible, 
then we can conclude that the current-voltage characteristic will fit the 
ideal formula. For these assumptions 8<p Kp on Fig. 5 will be small and the 
principal change in <p p will occur relatively deep in the w-region, at least 
beyond the transition region. So long as the hole concentration introduced 
in the w-region is much smaller than n„ , the hole current into the w-region 
will be a linear function of the value of p at the right edge of the transition 
region, being zero when p equals p n , the equilibrium value of p. This 
leads at once to a hole current proportional to p — pn and since the shift 
of <p p in respect to ^ at the edge of transition region is 5<p, p — pn is 
equal to p n (exp(a8<p/kT)-l). (These ideas are discussed in detail in Sec- 
tion 4.) A similar relationship will hold for electrons entering the ^-region; 
hence the total current will vary as exp(qd<p/kT)-l. This is a theoretical 
rectification formula" giving the maximum rectification for carriers 
of charge q. 

4. Treatment of Particular Models 

4.1 Introduction and Assumptions 

In this section we shall deal chiefly with good rectifiers so that the IR 
drop, discussed in connection with Ri in Section 3, is negligible. We shall 
deal chiefly with the case for which the transition region is narrow com- 
pared to the diffusion length; consequently, there is little change in I p in 
traversing the transition region. In Fig. 6(a) we consider a hypothetical 
junction in which the properties are uniform outside the transition region. 
The division of the specimen into three parts as shown is seen to be reason- 
able for germanium: In w-type germanium, the diffusion constant for holes 
is about 40 cm 2 /sec and the lifetime is greater than 10" 6 sec so that the 
diffusion distance is L v = VD^, > 6 X 10" 3 cm. This is much greater than 
most transition regions. 

The major drop in <p p must occur to the right of the transition region. This 
follows from our assumptions: First, we may neglect the IR drop in the 
/^-region; hence ip p is substantially constant from x = x a to * = x Tp . Second, 
the decrease in <p p is much less in the transition region than in the w-region; 
this follows from two considerations: the resistance for hole flow is lower in 

»C. Wagner, Phys. Zeits. 32, 641-645 (1931). 
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the transition region than in the w-region; the effective length of flow in the 
w-region, being L p , is greater than the width of the transition region. Con- 
sequently, the variation of ip p shown in Fig. 6(c) is seen to be reasonable. 
Similar considerations apply to <p n . As is shown in Fig. 6(c), the application 
of b(p does not alter ip p — \p in the /^-region nor <p n — ^ in the -w-region. The 
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Fig. 6 — Simplified model of a p-n junction. 

(a) Distribution of donors and acceptors. 

(b) Potentials for thermal equilibrium. 

(c) Effect of 5^ applied potential in forward direction. 



reason, as discussed in connection with (2.31), is that in these regions elec- 
trical neutrality requires an essentially constant value for the more abundant 
carrier. Hence the relationships between the tp's and ^ follow from (2.4). 

The nature of the potential distribution in the transition region has no 
effect in the considerations just discussed. However, as shown in Section 2, 
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the capacity of the transition region, which we shall denote by C T in this 
section, does depend on the nature of the transition region and, consequently, 
on the value of A". 

If the sizes of the ^-region and //-region are large compared to the diffusion 
lengths, we may assume the current at x a to be substantially I p only and 
that at x b , /„ only. The total current entering at x a can be accounted for 
as doing three things: (1) neutralizing the electron current flowing into the 
^-region across x Tp , (2) contributing to the charge in the transition region 
(this corresponds to the capacity discussed in Section 2) and (3) contributing 
a current flow to the right across x T „ ■ 

We have selected the hole current for analysis because the hole has a 
positive charge and the connection between the algebra and the physical 
picture is simplified. For the same reason, the text emphasizes forward 
current, although the equations are equally applicable to reverse currents. 
Nothing essential is left out by this process; since the sample as a whole 
remains uncharged, the current / is the same for all values of x and if I p is 
known, then I„ = I — I„is also determined. 

4.2 Solution for Hole Flow into the n-region 

We shall calculate first the hole current I p {x Tn ) flowing across x Tn . It is 
readily evaluated as follows: The value of p(x Tn ) is given by 



. aiflkT 

= p n e 



(4.1) 



where p n is the hole concentration in the n-region for thermal equilibrium. 
If we apply a small a-c. signal superimposed on a d-c. bias so that 

&p = v + w" (4.2) 

where vi is an a-c. signal, assumed so small that linear theory may be em- 
ployed (i.e. 7'i « kT/q), then 

P(xm) = (p n e qvolkt )(l + foft/WV). 

We resolve this density into a d-c. component po and an a-c. component p\ 



e 



p(xm) = pn + fr + pie ial (4.3) 
where 

h = p n {e^ ,kr -\) (4.4) 

Pi = (qP*v l /kT)e q ">" eT . (4.5) 

So long as p(x T „) « n n , the normal concentration of electrons in the 
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//-region, the lifetime r ,, and diffusion constant D for a hole will be sub- 
stantially unaltered by V- Application of the hole-current equation to the 
hole density p(x, I) gives 

I r =-qD d J.. (4.6) 

ox 

Combining this with the recombination equation 

| = t±Zl _ I «l = t±Z± + D £| (+ . 7) 

at t p q ox t p ax* 

leads to the solution 

P = p n + ^Tn^VDTp + ^ e i«W»: r .-)a+f«r 1 ,)l/»/(I>r J ,>l/» (4 g) 

The quantity \/^rp is the diffusion length and is denoted by Z p . (We shall 
use subscript /» for holes in the //-region and n for electrons in the ^-region 
for both L and t.) 

When p is large compared to p n , but small compared to n n , the ex- 
pression for p leads to the following formula for <p p : 

<P P =<Pn + v c - (kT/q)(x - Xm)/L p + i^-^tW+W-Mp, ( 4 . 9 ) 

This shows that the d-c. part of <p p varies linearly in the w-region, for large 
forward currents, and decreases by (kT/q) in each diffusion length L p . 
The transition from this linear dependence to an exponential decay for <p p 
comes when^p — tp n = (kT/q). This behavior of the d-c. part of <p p is useful 
in connection with diagrams of <p p versus distance. (See Sections 5 and 6.) 
The solution just obtained for p gives rise to a current at x Tn of 



ox 

= qp D/L p + qp,De tuT {\ + tor,)" 2 /!, 
The d-c. part is calculated by substituting (4.4) for p : 
I p0 (xm) = (qPnD/L P Ke' IVolkT - 1) ; 
= I p »{e q " olkT - 1) 
and the a-c. part is similarly obtained from (4.5) for pi : 

TAxm) =(qPnn/Lp)\e<" rn,kT) }(\ + iW,) 1 "* e M 
= (G p + iS p )vi6 mt = A p vie ut 



(4.10) 



(4.11) 



(4.12) 



where A p is called the admittance (per unit area) for holes diffusing into the 
//-region; its real and imaginary parts are the conductance and suscept- 
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ance. For wt p small, the real term G p is simply conductance per cm 2 
of a layer L p cm thick with hole conduction corresponding to the density 
p n + f ; it is also the differential conductance obtained by differentiating 
(4.11) in respect to !q . For the case of zero bias this establishes the result 
quoted in Section 1 that the voltage drop is due to hole flow in the w-region 
where the hole conductivity is low. 

In this section we hive treated r p as arising from body recombination. 
In a simple whose y and z dimensions are compirable to L p or L n , surface 
recombination may play a dominant role. However, as we show in Appendix 
V, the theory given here may still apply provided appropriate values for 
t p and t„ are used. 

4.3 D-C. Formulae 

The total direct hole current flowing in at x a is I p0 plus the current re- 
quired to recombine with electrons in the //-region. This latter current is, of 
course, equal to the electron current flowing into the //-region. This electron 
current, denoted by I,o or /,.o(*r P ), is obtained by the same procedure 
as that leading to (4.11) for I l0 except that bD replaces D and the subscripts 
of L and t are now n. Combining the two currents leads to the total direct 
current: 

+ / n0 = ( qC ) (fr 4- *&) (e" VolkT - 1) (4.13) 

for the direct current per unit area for applied potential v . The algebraic 
signs are such that I > corresponds to current from the //-region to the 
M-region in the specimen; r > corresponds to a plus potential applied 
to the //-end. The ratio of hole current to electron current across the transi- 
tion region is 

VbD7n 
VDt~ p 



h = I 



pO 




(4.14) 

_ Pf . /bjln = . Ap 
bl'n 'V Pp V '" 

where we have used the relationships n„p n = n p p p = fl,- from (2.2) and 
T p iin = T„p p = l/f from (3.2) and (3.3). These results can be summarized 
by saying that the current flows principally into the material of higher re- 

12 For convenience we repeat the definitions here: q = magnitude of electronic charge; 
D= diffusion constant for holes; p„ and n„ = thermal equiliurium value of p and «, as- 
sumed constant throughout //-region (.r > x T „); n,. and p p = similar values for .r < x Tp ; 
L p = dilTusion length = y/Ur p for holes in «-rcgion; r v = lifetime of hole in M-region be- 
fore recombination; b = electron mobility/hole mobility; L„ and t„ similar in quantities 
for electrons in /"-region; a n = qnbn,, and a p = qnp p are the conductivities of the two regions. 
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sistivity. We can also say that the hole current depends only on the >?-type 
material and vice versa. For a p-n junction emitter in a transistor with an 
w-type base, it is thus advantageous to use high conductivity p-type ma- 
terial so as to suppress an unwanted electron current. 

For comparison with experiment, it is advantageous to express the values 
of p n and ?i p in terms of the conductivities a n and <r p . If the conductivity of 
the intrinsic material is written as 

<n = qimtil + b), (4.15) 

then, if p n <<C n n and n p <5C p p , we find 

qnp n = bal/(l + bfv n (4.16) 

qnbn p = b<xV(\ + bf<j p . (4.17) 

Using these equations, we may rewrite (4.11) and a corresponding equation 
for electron current into the /»-region so as to express their dependence ond-c. 
bias v Q and the properties of the regions: 



-/pc(flu) = 



Ino(vo) = 



ba'i kT , QVo lkT 


1) 


(i + b)*a n L p q kC 


m G p0 - ie"" '" 7 - 1) 




= I pa (e qv °' kT - 1) 




bOi kT , nvQlkT 


1) 


{\ + bYo p L n q K 


= G„o — (.e qv °' kT - 1) 
1 





(4.18) 



(4.19) 



= In,(e qvo ' kT - 1). 



The values of G p o and Co (which are readily seen to be the values of the low- 
frequency, low-voltage (ro < kT/q) conductances) and the saturation 
reverse currents are given by 

G *» s JXTT)^T P s kf Ipl {m 

G "° s (TTW^Tn s It Ina (4 - 21) 

The expression for direct current then becomes 



/•GO - [^ + Gj(-)[e" 0/lr - 1] 

= (/„. + in.)[<r o,kT - i]. 



(4.22) 
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4.4 Tolal Admittance 

In order to calculate the alternating current, we must include the capacity 
of the transition region, discussed in Section 2. Denoting this by C T , we 
then find for the total alternating current. 

lac = (G p + iS p + G n + iS n + faCr) Pi = Avx (4.23) 

where G n and S n are similar to G p and S p but apply to electron current into 
the /»-region. The value of the hole and electron admittances can be ex- 
pressed as 

A v = G p + iS p - (1 + iuT p ) U2 G f oe 9V o lkT (4.24) 

A n ■= G n + tf„ = (1 + io 3Tn ) in G nQ e <lVolkT (4.25) 

For low frequencies, such that u> is much less then 1/r, , we can expand 
G p + iS p as follows: 

G„ + iS P = G+r*" + utrJliG+r** (4.26) 

Hence (r p /2)G j0 c 8 * '*' behaves like a capacity. 

It is instructive to interpret this capacity for the case of zero bias, t-o = 0, 
for which we find: 

C p = r p G p0 /2 = T p gp n n/2L P = c } 2 p n L p /2kT. (4.27) 

The last formula, obtained by noting that t p /j. = qrpD/kT = gL p /kT, has 
a simple interpretation: gp n L p is the total charge of holes in a layer L p thick. 
For a small change in voltage v, this density should change by a fraction 
qv/kT so that the change in charge divided by the change in v is 
(q/kT)(qp n L p ) which differs from C p only by a factor of 2, which arises from 
the nature of the diffusion equation. 

This capacity can be compared with C T nm%. , discussed in Section 2, (see 
equation (2.39) and text for (2.42)) for germanium at room temperature 
as follows: 



C p _ q~ p n L p kTa _ pnL p a 
CV^t. ~ 2kT ' io 5 2 n\ 20nf " 



(4.28) 



For a structure like Fig. 6(c), the excess of donors over acceptors reaches its 
maximum value, equal to />„ , at x Tn leading to «„ = ax Tn . Consequently 
a = njxrn ■ Substituting this value for a in (4.28) and noting that 
pnii„ = in gives 



Cp _ L 



Ct nffllt. 20.T7 



(4.29) 



As discussed at the beginning of this section, L p — 6 X 10 cm for holes 
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in germanium. Hence if the transition region is 6 X 10~ cm thick, the 

diffusion capacity C ,, will dominate the capacitative term in the admittance. 

Although A p simulates a conductance and capacitance in parallel at low 

frequencies, its high-frequency behav'or is quite different. In Fig. 7 the 
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Fig. 7 — Real, G, and imaginary. S, components of admittance for hole flow into w-region. 

(a) MPAp/Gpo = 10 3 (1 + /ur,,) 1 ' 2 corresponding to uniform w-region. 

(b) 10 2 X Formula of Appendix III. corresponding to layer of high recombination 
rate in front of w-region. This causes G to exceed S at higher frequencies than for (a). 

(c) 10 X Equation (4.33), corresponding to a retarding held in the w-region, with 

L T = Lp/VlO- 

(d) Equation (4.33) with U = L P /10- 

behavior of (1 + ia)T p ) 112 = A v /G r ,o , is shown. For high frequencies G p and 
S p are equal: 



Gp — S p — v Tp/2 G P o vto — 



(1 + bfan V2D 



(4.30) 



Thus for high frequencies the admittance is independent of r„ and is deter- 
mined by the diffusion of holes in and out of the H-region. The three straight 
asymptotes have a common intersection at the point G p o , or = 2 on Fig. 7, 
a fact which is useful in estimating the value of r from such data. 

For large o>, S p varies as u> as shown in (4.30) whereas St is wCr • Hence 
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at very high frequencies C T will dominate the admittance. At very high fre- 
quencies C T itself will have a frequency dependence; however, for the as- 
sumptions on which the treatment of this section is based, the relaxation 
time for the transition region t t is much less than t p . This is a consequence 
of the fact that, although diffusion of holes into the transition region is 
required for the charging of C T , the distance is relatively short, being in 
fact only that fraction of the width x Tn — x Tp of the transition region in 
which \f/ rises by kT/q; in germanium this will be about one-tenth of 
*m — Xt p • Since diffusion times vary as (distance) 2 , the ratio of the times is 

Tt ._ \Xtti — X pn ) (431) 

t p lOOLp 

Hence if L p > x Tn — x pn , t t will be much less than r p . 

4.5 Admittance Due to Hole Flow in a Retarding Field 

In Appendix II we treat the case in which a potential gradient, due 
to changing concentration for example, is present in the n- and ^-regions. 
This tends to prevent holes from diffusing deep in the w-region and for 
this reason the w-region acts partly like a storage tank for holes under a-c. 
conditions, thus enhancing S p compared toG p in A p . If the electric field is 
—d\J//dx = kT/qLr , where L r is the distance required for an increase of kT/q 
of potential (i.e. a factor of e increase in »„), then the value of A p is 

A p = [q» Pn /L p ] —^—j— . WTp)(2jLr/£p)2]1/2 l«4 

For ut p > 1, this admittance is largely reactive provided 2L T /L p is suffi- 
ciently small. 

The dependence of A p upon co is shown in Fig. 7 for two values of L r /L p . 
The plot shows the real and imaginary parts of 

A p /[2qnp n Lr/L 2 p ] = 1+[1 + (1 + iwT J (2£f/Zj))2 ,i, s ( 4 - 33 ) 

for L p /Lr = 10 1/2 and L p /L r = 10, the two curves being relatively displaced 
vertically by one decade. The second value implies that the field keeps the 
holes back so that they penetrate only xfr their possible diffusion length in 
no field. It is seen that for this case the storage effect is very pronounced and 
the susceptance S is much larger than G for high frequencies. 

The function (1 + k»>t„) 1/2 , discussed earlier, corresponds to the limiting 
case of (4.32) for L T = * . 

13 In Appendix IV an analytic treatment of Ct is given. 
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4.6 The Effect of a Region of High Rale of Generation 

There is evidence that imperfections, such as surfaces and cracks, add 
materially to the rate of generation and recombination of holes and elec- 
trons. If there is a localized region of high recombination rate in the transi- 
tion region, there will be a pronounced modification of the admittance char- 
acteristics. In Fig. 8(a) such a layer is represented at * = 0. In Fig. 8(b) 
the customary plot of <p p and tp„ versus x is shown. If we neglect the effect 
of the series resistance terms denoted by Ri in Section 3, the change Sip will 
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Fig. 8 — The effect of a localized layer of high recombination rate on the junction 

characteristic. 

(a) Location of layer of high recombination rate. 

(b) Quasi Fermi levels. 

(c) Distribution of hole current showing rapid change at layer of high recombination 
rate. 



occur in the ^-region for <p„ and in the w-region for <p p . The hole current flow- 
ing into the //-region will thus be the same as before and will be given by 
equation (4.11) or (4.18) and denoted by Ipo{Sip). Similarly, the electron 
current will be I n o{Sip). In the layer we shall suppose that there is a rate of 
generation of hole electron pairs equal to g a per unit area of the layer and a 
rate of recombination proportional to r a np per unit area. We suppose, further- 
more, that the layer is so thin that // and p are uniform throughout the layer. 
The net rate of generation is thus 



g a - r a np = g a [l - e ql *>-^ ,kT \ 



(4.34) 
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since for equilibrium conditions the rates balance so that r a n\ = g a . The 
net hole current recombining in the layer per unit area is thus 

Ir(* P -<Pn) = qg a [e^"-^ lkT - Y] (4.35) 

There must, therefore, be a discontinuous decrease of hole current across 
the layer. The total hole current flowing in at x = x a , which is also the total 
current /, thus does three things: for x < x Tp , it combines with I„ (8<p); 
for x Tp < x < xT n , it combines with electrons at rate I r (&p)l for % > x T „ , 
it flows into the «-region in amount Ipo(fo). This leads to 

In other words the layer of high recombination acts like a rectifier 
in parallel with 7 n0 (M + /*o(M- The frequency characteristic of I r (ty), 
however, will be independent of frequency and will contribute a pure con- 
ductance to the admittance of the junction. 

If the layer is considered to have finite width, however, it will exhibit 
frequency effects just as does I p in the //-region. In Appendix III, we treat 
a case in which the layer is a part of the //-region itself but has a recombina- 
tion time different from the main layer. If the time is shorter, a large amount 
of the hole current may recombine in this layer. For high frequencies, the 
current may not penetrate the layer, in which case the admittance for hole 
current is determined by the thin layer rather than by the whole ;/-type 
region. A case of this sort is shown in Fig. 7. In this case the thickness of the 
layer is \ of its diffusion length and in it the lifetime of a hole n is £ the 
value t p in the main body of the //-region. The hole current will thus be 
restricted to this layer when the diffusion distance -y/D/u is less than the 
layer thickness (£) \/Dti; this corresponds to wr t > 9 or ut p > 81. The 
presence of the high rate of combination in the layer is evidenced by the 
tendency of G to be greater than S at high frequencies. If the layer were 
infinitely thin, as discussed above, it would simply add a constant conduct- 
ance to the admittance. 

4.7 Patch Effect in p-n Junctions 

If there are localized regions of high recombination rate, a "patch effect" 
may be produced in an n-p junction. As an extreme example, suppose the 
value of g a for the layer just considered is allowed to become very large; then 
the recombination resistance may become small compared to Ri in Section 3 
and the junction will become substantially ohmic. If the region of high 
rate of recombination is relatively small compared to the area of the rest 
of the junction, then the behavior of the junction as a whole may be re- 
garded as being that due to two junctions in parallel. Over most of the area, 
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the currents will flow as if the patch were not present so that one compo- 
nent of the current will be that due to the uniform junction. In addition 
there will be current due to recombination and generation in the patch. 
The series resistance to the patch will be relatively high due to the constric- 
tion of the current paths. On the other hand, the value of 7 r (50) associated 
with the patch may be very high. Hence the current due to the patch will 
be that of a low impedance ideal rectifier in series with a high resistance; 
and if the ratio of impedances is high enough, such a series combination 
amounts essentially to an ohmic leakage path. Thus patches in the p-n 
junction will tend to introduce leakage paths and destroy saturation in the 
reverse direction. 

An extreme example of a region of high rate of recombination would be a 
particle of metal making a non-rectifying contact to both p- and w-type 
germanium. Since holes and electrons are essentially instantly combined in a 
metal, the boundary condition at the metal surface would be equality of 
<p p and (p„ . This would mean that near the metal particle, p p and <r„ could 
not differ by 5^, the condition required, over some parts of the junction at 
least, in order for ideal rectification to occur. 

A common source of imperfection in p-n junctions arises from dirt or 
fragments on the surface which overlap the junction. Even if these do not 
actually constitute a short circuit across the junction, they may furnish 
patches of the sort discussed here and modify the junction characteristic. 

4.8 Final Comments 

Another possible cause for frequency effects may be found in the trapping 
of holes or electrons." When an added hole concentration is introduced into 
an w-region, a certain fraction of the holes will be captured by acceptors and 
later re-emitted or else recombined with electrons while trapped. Investiga- 
tion of this process is given in Appendix VI. One interesting result is that the 
trapping of holes in a uniform »-region cannot produce an effective suscep- 
tance (i.e. iwC) in excess of the conductance, as can a retarding field. 

Finally it should be remarked that important and significant variations 
of the conductivity in the p- and w-regions may be produced by hole or 
electron injection. Under these conditions, when the hole concentration 
approaches n„ ,\f/ — <p n will vary. Under these conditions Ri may be appreci- 
ably altered. These factors favor the p-n junction as a rectifier since they lead 
to a reduction of series resistance under conditions of forward bias and thus 
tend to improve the rectification ratio. 

14 Frequency dependent effects in C112O rectifiers have been explained in this way by 
J. Bardeen and W. H. Brattain, personal communication. 
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5. Internal Contact Potentials 

The theory of p-n junctions presented above has interesting consequences 
when applied to the distribution of potential between two semiconductors 
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Fig. 9— Internal contact potentials showing how presence of injected holes produces a 
contact potential across /». 

under conditions of hole or electron injection. In Fig. 9 we illustrate an 
X-shaped structure. A forward current flows across the junction Pi and 
out of branch Ni . If the distance across the intersection is comparable with 
or small compared to the diffusion length for holes, a potential difference 
should be measured between P 2 and N 2 . The reason for this is that holes 
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flow easily into P% since the potential distribution there favors their en- 
trance. Since, however, P 2 is open-circuited this hole flow biases J« in the 
forward direction; since J* is high resistance, an appreciable bias is developed 
before the counter current equals the inward hole flow and a steady state is 
reached. No similar effect occurs in the branch A r 2 ; consequently Pi will be 
found to be floating (open-circuited) at a more positive potential than A r 2 . 

Parts (b) to (e) describe this reasoning in more complete terms. We 
suppose that the ^-regions are more highly conducting than the ^-regions 
so that the current across J\ , shown in (b), is mainly holes. The potentials 
(f p and (f n along the .r-axis will be similar to those of Figs. 5 and 6; (c) shows 
this situation and indicates that the diffusion length for electrons in the 
/»-region is less than for holes in the M-region. Along the y axis <p p and <p n 
vary as shown in (e), the reasoning being as follows: At the origin of coordi- 
nates <p p and ip n have the same values as for (c). The transverse hole current 
(d) has a small positive component at y = since, as mentioned above, P 2 
tends to absorb holes and thus increase diffusion along the plus y-axis. Since 
the net transverse current is zero, /„ = —I p in (d). The <p curves of (e) 
have been drawn to conform to the currents in (d) ; <p„ is nearly constant 
in the M-region and <p p is nearly constant in the ^-region. As concluded in 
connection with Figs. 5 and 6, <p n and ip p are also nearly constant across the 
transition region. These conclusions lead to the shape of <p n and <p p for y > 
in (e). For y < 0, the reasoning is the same as that used in Sections 3 and 4 
and we conclude that <p n is essentially constant. Hence, a difference in the 
Fermi levels at Po and N% will result. 

In Fig. 10 we show a structure for which we can make quantitative calcu- 
lations of the variations of ip p and <p n . We assume for this case that the 
forward current from Pi to N does not produce an appreciable voltage drop, 
i.e. change in \p and <p„ , in region N. This will be a good approximation if the 
dimensions are suitably proportioned. We shall next solve for the steady- 
state distribution of p subject to the indicated boundary conditions assuming 
that p is a function of x only. As we have discussed in Section 4.1, when p is 
small compared to n in the M-region, we can write 

p = p n e ^p-.„mr (5>1) 

In keeping with the treatment in the next section of this structure as a 
transistor, the terminals are designated emitter, collector and base, the po- 
tentials with respect to the base being tp e and <p c . The contact to N or the 
base is such that tpt ■= <p n in this region. Hence, the boundary conditions at 
J\ and J 2 are 

p x = p n e q * JkT x = -w (5.2) 

p* = p„e q *' lkT x=+w (5.3) 
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The function p{x) which satisfies these boundary conditions and the equation 

D £p _ P-Pn m q (54) 

dx 2 t p 

is 

FW ^ 2 cosh (w/L p ) 2 sinh (w/Zy) 

(5.5) 




v = V- + v, 



v = v c +*> b 



■x a -w o w 

DISTANCE THROUGH SAMPLE X— *■ 

Fig. 10— Model used for calculation of internal contact potential and to illustrate p-n-p 

transistor. 

(a) Semiconductor with two p-n junctions and ohmic metal contacts. 

(b) Quasi Fermi levels showing internal contact potential between b and c. 

which gives rise to a hole current across J* into Pi of amount 

dp I 

/p= ~ qD £\ 

QX \x=w 

= q £~ \(pi - P*) coth f + (2/> n - * - fc) tanh f 1 
2L P L ^p -^p- 1 

- (/>-> - />») (coth ~ + tanh £"- j 

ft,^r [e w/ * r - 1] _ (e w ' /tr - 1) 1 
L v Lsinh (2w/L p ) tanh {2w/L p )\ 

= cscli (2w/L p )Ipo(<p t ) - coth (2wLp)I p o(<Pc) 



(5.6) 
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where, by I^o (<p), we mean the hole current which would flow in the forward 
direction across either Ji or Ji if uninfluenced by the other (i.e. the function 
of (4.11) or (4.18) and (4.20).) The equation shows that a fraction csch 
(2w/L p ) of the current 7 p n (<p t ), which would be injected by <p f on i\ in the 
absence of J-> , flows into P* . The conductance of P 2 across 7 2 is increased 
by the factor coth(2w/Z„). 

The current into Pi carried by electrons will be unaffected by J\ and can 
be denoted by — I„o(<p c ) the minus sign resulting from the fact that currents 
into P-i are in the reverse direction. The total current flowing into P 2 contains 
the —I n o(<p c ) and —I po (<f>c) terms and must cancel the -\-I p o(<p e ) term for 
equilibrium. Hence: 

/ n o(vc) + coth (2w/L p ) I p0 (<p c ) = csch (2w/L p ) Jpofo) (5.7) 

If p n )j> n p , the I„o term can be neglected compared to coth (2w/L p ) I p0 . 
Hence the value of <p c must satisfy 

I P o(<Pc) = sech (2w/L p ) I p0 (<p,). (5.8) 

For <p f > kT/g, the exponential approximation may be used for I p o in both 
terms: 

<p c = * - (kT/q) In cosh (2w/L p ), (5.9) 

so that, if (2w/Lp) is the order of unity, ip c should be only about (kT/g) 
less than p t . For (2w/L p ) large, we get 

& = <p. - (kT/g) (2w/L p ) (5.10) 

corresponding to the linear drop of <p p , discussed in connection with equation 
(4.9), across the distance 2w. 

When if e is negative, so that we have to deal with reverse current, p c 
will not decrease indefinitely but will reach a minimum value given by 

[espgjp,/kT\ - 1 = -sech (2w/L p ) (5.11) 

and corresponding to saturation reverse current across J\ , so that 

<p e = -(kT/g) In [1 + (1/2) csch*(w/L p )]. (5.12) 

The floating potentials of p-type contacts to »-type material into which 
holes have been injected (or w-type contacts to p-type material with in- 
jected electrons) are reminiscent of probes in gas discharges which tend to 
become charged negative in respect to the space around them because they 
catch electrons more easily than positive ions. The situation may also be 
compared with that producing thermal e.m.f.'s; in fact a "concentration 
temperature" of the semiconductor with injected holes can be defined by 
finding the temperature for which np = >ii(T). We conclude that, in the 



472 BELL SYSTEM TECHNICAL JOURNAL 

absence of thermal equilibrium, different potentials depending on the nature 
of the contact are, in general, the rule rather than the exception. 

The bias developed on P 2 or c will change its conductance. If we suppose 
that <p t and <p b are held constant, then the current flowing into c is obtained 
by the same reasoning that led to (5.7) and is 

Ic(<Pc, <f t ) = iAfc) + coth — ipoOc) - csch -- I P oM- (5.13) 

L p Lp 

For an infinitesimal change in <p c from the value which makes I c {<Pc , <p«) 
vanish, the admittance to c is readily found from (4.18) and (4.19) to be 

( -— ° ) = l'no(<Pc) + COth j- l' p o(<Pc) 

Kd(pc/ *' (5.14) 



= GnO + COth -- GpO 



e QV c lkT 



which shows that pronounced variations in admittance should be associated 
with variations in hole density in N in Fig. 10. 

6. p-n-p Transistors 

The structure shown in Fig. 10 is a transistor with power gain provided 
the distance w is not too great. As a first approximation, we shall neglect 
the drop due to currents in the N region. If we use P 2 as the collector and call 
the collector current, L , positive when it flows into P 2 from outside, we shall 
have from (5.13) 

I e = -csch y^ Ipoift) + coth -^ I v n(<Pc) + /»o(*»e). (6.1) 

Lp L p 

The emitter current is similarly 

I, = coth j- IgoM - csch — IpoM + ImM- (6.2) 

L p L p 

If p n » n p , then the 7„o terms can be neglected. However, the base current 
will not vanish but will be 

Ib = -I ( - I c = [csch ~ - coth ^ 1 [IM + iAfc)} 

1 L " LvJ (6.3) 

= 2sinhW£, { ^ + lMc)l 
sinh Iw/Lp 

16 The variations in admittance discussed in connection with metal point contacts in 
an accompanying paper in this issue (W. Shockley, G. L. Pearson and J. R. Haynes, Bell 
Sys. Tech. JL, July, 1949), arise from this cause; however, the nature of the contact is 
not as simple as here. 
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For w/Lp large, the junctions do not interact and the hyperbolic coefficient 
becomes unity and I b = —[T p0 (<p f ) + /»i(ft)]. 
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DISTANCE THROUGH SAMPLE , X- 



(a) Thermal equilibrium. 

(b) Operating condition. 



Fig. 11 — p-n-p transistor. 



If <p c is several volts negative, so that I p o(f c ) has its saturation value 



I „ a (see (4.11) and (4.20)), then the ratio —8I r /8T t = a has the value 



u. csch r, . 2w 
a =-sr. = — ^ = sech z;- 

coth j- 

J-'-n 



(6.4) 



For (2w/L p ) = 0.5, 1, 2 respectively, a = 0.89, 0.65, 0.27. Since the output 
impedance R>>o will be very high when <p c is in the reverse direction, and the 
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input impedance will be low, the power gain formula 6 a R-a/Ru will yield 
power gain even when a is less than unity. 

In certain ways the structure of Fig. 10 resembles a vacuum tube. In Fig. 
11, we show the energy band diagram, with energies of holes plotted upwards 
so as to be in accord with the convention for voltages, (a) shows the thermal 
equilibrium distribution and (b) the distribution under operating condi- 
tions. It is seen that the potential hill, which holes must climb in reaching 
the collector, has been reduced by <p t . The n-region represents in a sense the 
grid region in a vacuum tube, in which the potential and hence plate current, 
is varied by the charge on the grid wires. Here the potential in the //-region 
is varied by the voltage applied between base and emitter. In both cases one 
current is controlled by another. In the vacuum tube the current which 
charges the grid wires controls the space current. Because the grid is negative 
to the cathode, the electrons involved in the space current are kept away from 
the grid while at the same time the electrons in the grid are kept out of the 
space by the work function of the grid (provided that the grid does not 
become overheated.) In Fig. 11, the electrons flowing into the base control 
the hole current from emitter to collector. In this case the controlled and 
controlling currents flow in the same space but in different directions because 
of the opposite signs of their charges. 

As this discussion suggests, it may be advantageous to operate the p-n-p 
transistor like a grounded cathode vacuum tube, with the emitter grounded 
and the input applied to the base. 

The p-n-p transistor has the interesting feature of being calculable to a 
high degree. One can consider such questions as the relative ratios of width 
to length of the u-region and the effect of altering impurity contents and 
scaling the structure to operate in different frequency ranges. However, 
we shall not pursue these questions of possible applications further here. 

Acknowledgment 

The writer is indebted to a number of his colleagues for stimulating dis- 
cussions and encouragement, in particular to H. R. Moore, G. L. Pearson and 
M. Sparks, whose experimental work, to be described in later publications, 
suggested development of the theory along the lines presented above. He is 
also indebted to J. Bardeen, P. Debye, G. Wannier and W. van Roosbroeck 
for theoretical comments and suggestions and especially to the last and 
to Mrs. G. V. Smith for valuable assistance in preparing the manuscript. 

10 Physical Principles Involved in Transistor Action, J. Bardeen and W. H. Brattain, 
Phys. Rev. 75, 1208 (1949). 



p-n JUNCTIONS IN SEMICONDUCTORS 475 

APPENDIX I 

A Theorem on Junction Resistance 

We shall here prove that the junction resistance is never less than the value 
obtained by integrating the local resistivity l/qnip + bn). This is accom- 
plished by analyzing the following equation which we shall discuss before 
giving the derivation: 

the meaning of the symbols being that shown in Fig. 5. This expression is 
valid even if large disturbances in p and n from their equilibrium values 
occur. The second integral is positive since the integrand is never negative. 
It may be very large if <p p — <p n ^> kT/q in some regions. If, in the first in- 
tegral, we consider that I p and /„ may be varied subject to the restraint 
Ip + In = I, we may readily prove that the first integrand takes on a mini- 
mum value when 

pi bnl 

Ip = —— — — and I n = 



p + bn p -f- bn' 

For this minimum condition, the first integral becomes 

f I dx/qn{p + bn) = J*R> 

where Ro is simply the integrated local resistivity. If / does divide in this 
way, the second integral is zero, a result which we can see as follows: 

/, = —qupdipp/dx 

I„ = —qpbn dtp„/dx 

d<p p /dx _ Ip/p 
d(p„/dx In/bn 

Hence, if the current divides in the ratio of p to bn, then d<p p = d<p n and, 
since ip p = <p n at x a , <p p = <p n everywhere and the second integral vanishes. 

In general, of course, the conditions governing recombination prevent 
current division in the ratio p'.bn and then Sip /I > Fq • 

The equation discussed above is derived as follows: We suppose that 

Pj»(#a) = P»(#a) = <Pa 
<P P (*b) = <Pn(-\'b) = <Pb " 
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Then 

18<p = — (Ip <fp + In <Pn) 

x a 

f Zb d 
J Xa dx 

-jC(S*+S-)*-jC(^+^)* 



Since 



and 



dx dx 



these two integrals are readily transformed into the ones previously dis- 
cussed. 

APPENDIX II 

.Admittance in a Retarding Field 

We shall here derive the admittance equation for holes diffusing into a 
retarding potential \p = kTx/qL r in which the potential increases by kT in 
each distance L r . The differential equation for the a-c. component of p is 

^--r r -^[- D d -*=J- 

This equation may be solved by letting p = pi exp(iat - 7*) as may be 
seen by rewriting the equation and substituting this expression for p: 

Bf^ + fafl-io + w* 

= - 7 Z? J -7 + i] £ - i (1 + *wrp)# « 

leading to 

1 + [1 + (2L r /L p )\l + i«r p )] 1/2 



7 = 



2L r 
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The corresponding current evaluated at x = where p = pi exp(iwt) = 
(p „qvi/ kT)exp(iwt) is given by 

/--«r»3?+*a 



a.r 



^(1 + iwT p )/> 



» + i 1 ; 



_ p»q\l + tut) 2L r 



Vie 



*7Y P 1 + [1 + (2L r /L p y(\ + tor,)] 1 ' 2 

_ q np n 2L r (1 + io>T p ) M 

~ L\ " 1 + [1 + (2L r /L p ) 2 (l + *W P )] ,/2 ' Vl€ 

= A p vie tul . 

This is equivalent to (4.32) in Section 4. 

APPENDIX III 

Admittance for Two Layers 

We shall here treat a case in which there is a thin layer on the n-side of 
the transition region in which recombination occurs much more readily 
than deeper in the w-layer. The case of an infinitely thin plane, discussed in 
Section 4, is a limiting case of this model. We shall suppose that the layer 
extends from x = — c to x = while x > corresponds to the w-region. We 
shall suppose that the potential in the layer is uniform with value \}/i whereas 
in the w-region it has value ^2 • The lifetimes of holes will be taken n and t 2 
in the two layers. The solutions for pi and pi are evidently 

pi = pio + (4 e~ ax + B e +ax ) e iat x < 



P2 = p*> + C e~**** x > 



where 



a = (1 + «#n) w /VKl - (1 + U*rd W /Li 

- (1 + W7V5S « (1 + ^r 2 ) U2 /L 2 . 

The boundary condition for continuity of <p p , required to avoid singularity 
in d<p p /dx, is 

p i e 9 * l,kT = p x e q * llkT 
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and, for continuity of hole current, is dpi/dx = dpt/dx. Expressing these 
in terms of A, B, C, a and for the a-c. components yields: 

A + B = Ce gl * l ~*' ),kT m CF 

a(A - B) = $C 

so that 

A = (F + (3/a)C/2. 

B = (F - f3/a)C/2. 

Hence the ratio — [dp/dx]/p at x = —c is 

d In p _ a(Ae +ac - Be~ ac ) = a(Fa sinh ac + g cosh ac) 
~~ ~~dx (A e +ac + Be~ ac ) Fa cosh ac + j3 sinh ac 

Since at » = — c, the a-c. component of pi is {qvi/kT)pvfi , the admittance is 

= -gg^ = (fDp u /kTK-d]np/dx) 
Vie xat 

, / r \/< . • ^l/2 Fa sinh ac + f9 cosh ac 

= ( ^ lo/Ll)(1 + toTl) Fa cosh ac + ^ sinh ac' 

For c — * 0, this transforms into 

(qnpio/Li) (1 + /con) 1/2 /3/Fa = (qn(pio/F)/L s )(l + icor 2 ) 1/2 

which agrees with Section 4, since pw/Fthen corresponds to p n . 

If c/Li and F are not large, an appreciable amount of recombination 
takes place for .r > for low frequencies. Dispersive effects will then occur 
corresponding to t 2 . The a-c. will not penetrate to x = 0, however, 
if c(co/D) 1/2 » 1 and the dispersive effects will then be determined by n . 

The frequency-dependent part of the admittance, 

... Fa sinh ac + ft cosh ac 
(1 + ten) pa cQsh ^ + ^ ginh ac , 

has been computed and is shown in Fig. 7 for t p = r 2 , F = 1, n = t 2 /9 and 
c /£j = |. For these values about half the hole current reaches * = for 
low frequencies. As the time constant for diffusion through the layer is t„/81, 
as discussed in Section 4.6, the layer will act as a largely frequency-inde- 
pendent admittance well above the point for wt p = 1. This is reflected in 
the behavior of the curves of Fig. 7 and, for frequencies in the \/ui range, 
it is seen that G is larger than 5 by about 50% of the low-frequency value of 
G; this split of G + iS into (£)G plus approximately (i)Go (1 + ia)T p ) m corre- 
sponds to the fact that about half the holes are absorbed in layer 1 for the 
assumed conditions. 
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APPENDIX IV 

Time Constant for the Capacity of the Transition Region 

For this case we shall consider the case of holes in an a-c. field with po- 
tential 



a V£/x 



(iwt\ 



where the d-c. retarding field is kT/qL r and the a-c. field is kT/qLi where 
l/Li is considered small for the linear theory presented here. The expression 
for the current of holes is 



dx OX 



-B+'& + 9] 



We shall obtain a solution for p by letting 

p — poe + pi [e — e \e , 

while neglecting recombination in this region so that p must satisfy the con- 
dition p = — d (hole current)/d.v leading to the differential equation 

*[&+2(£+0|->- 

There are three separate exponential dependencies of the variables leading 
to three equations (neglecting terms of order (l//.i) ) 

-xlLr+ivt. D K _ _ ^ _ _ p(> _ iwpi = Q 

e-" +iut : DW - y/L r ]px - u* - 

The first equation is satisfied by the equilibrium distribution and the 
second by 

pi = —po D/iu LiL r 

and the last by 

= 1 + Vl + 4'utf /D 
7 " 2L r 

It is evident that dispersive effects set in when 

co = D/4LI 
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This corresponds to the result used in (4.31) in which (a>„ - .v, P )/10 was 
used for L r . For smaller values of w the current may be calculated and put 
in simple form by expanding 7 up to terms including o> . The resulting ex- 
pression for the current is 

/ = —iu)qp L r (L r /Li)e l " 

This is interpreted as follows: The a-c. voltage across a layer L T thick is 

ty> = (kT/g) (Lr/L^ 1 

and, if we consider plus voltage as producing a field from left to right, then 
the a-c. voltage across L r is V = -5^. Substituting this for (L r /Li)exp(iut) 
gives 

/ = iaqp L r (q/kT)V 

Here qp L r is the total charge in the layer L T , (qV/kT) is an average frac- 
tional change in this charge for V so that (qp L r ) (qV/kT) -^ V is a capacity. 

APPENDIX V 

The Effect of Surface Recombination 

In this appendix we shall consider the effect of surface recombination upon 
the characteristics of the p-n junction. As for Section 4 we shall illustrate 
the theory for the case of holes diffusing into w-type material. For sim- 
plicity we shall treat a square cross-section bounded by y = ±w, z = ±w, 
the current flow being along -\-x. 

We shall denote the a-c. component of p as 

pi = px (*, y, z, t) 
At x = 0, the edge of the w-region, we shall suppose that <p p and \p are inde- 
pendent of y and z so that we shall have 

pi(0,y,*,t) = Pioe iut = (p n qoi/kT)e M 
by reasoning similar to that used for equation (4.5). The boundary condi- 
tion at the surface will be 

—D — = spi for y = -\-w 

dy 

This states that the recombination per unit area is spi and is equal to the 
diffusion to the surface —Ddpi/dy. Similar boundary conditions hold for the 
other surfaces. By standard procedures involving separation of variables 
we may verify that the solution satisfying the boundary conditions is 

CO 

pi = X G^e****** cos fay cos foz 

i.J-0 
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where the eigenvalues fr are determined by the boundary condition 

|8,ici tan /3,-w = sb/D = x- 

We use 0,- = (S{w for brevity later. Because of the symmetry of the boundary 
conditions it is not necessary to include sine functions in the sum. The value 
of an is given by 

an = (1 + iuTjtf/iDTjtf 

where r,/ is the lifetime of a hole in the eigenf unction cos /3,-y cosftz; i.e. 
th is the lifetime which makes 

p = exp (-l/tij) cos fry cos ft z, 

a function which satisfies the surface boundary conditions, a solution of the 
equation 

dp/dl = DV 2 p - P/t = -D($\ + $)p - pi ' T 

where to simplify the subsequent expressions we have omitted the subscript 
p from r. This equation leads to 

i - D{fi + fi) + -• 

Tif T 

The coefficients o,/ are readily found since the cos /3,y functions form an 
orthogonal set (as may be verified by integrating by parts and using the 
boundary conditions). The values are 

aa/pio = 4[sin 0,- sin 0,]/M,U + (1/2*0 sin 20,] -[1 + (1/20/) sin 20,] 

The current corresponding to this solution is 



h= ~qD f [ (dp/dx) dy dz 



integrated over the cross section at x = 0. This gives 

h = qDpi e lat 2 <Xii{a i j/pio)('hv/6idj) sin 0,- sin 0, 

Substituting for a,-,- and inserting />io = p n qi<i/kT, we obtain an expression 
for the admittance A p = Ii/Vi expO'w/): 

4 sin 2 0,- siii 2 0, 



-4 P = 4w qnpn 2,; a,-,- 



"'^[ 1 + fe) sin2 "'][ 1 + (i) sin29 '] 



where the sum plays the role formerly taken by (1 4- mty^/y/Dr in equation 
(4.12); the factor 4?c is the area of the junction. 

We shall analyze the formula for the case in which recombination on the 
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surface is smaller than diffusion to the surface so that x is not large. The values 
of Bi, over which the sum is to be taken, may be estimated as follows: in 
each interval of 0, of the form mr to (n + {$))*, 0. tan 0,- varies from to °c, 
giving one solution to t - tan 0, = x- For x small, the solutions are approxi- 
mately 

0o — sin O — tan O — Vx 

0i - 7r + xA; — sin 0i = tan 0i = xA 

n = mr + x/»7f; (—1)" sin 0„ = tan 0„ = x/ntr 
From this we see that the terms in the sum are as follows: 

«oo-4xVx 2 -i = «oo 
a n0 -2(x/nT) 2 /(mr) 2 = a n0 2x/mr 

,4/448 

From this it is evident that unless x is large, the series converges very 
rapidly. (This conclusion is not altered when the increase in <x„ m with /?„/?,„ is 
considered.) Thus the dominant term in the admittance is 

4w 2 qiipo (1 + hotoo) i/2 / \/Dtoo 
where 



1/roo = 2 (£) (0 2 o ) + 1/t 



^ 2 (?) ™ 

= 2 + 1/ 



: " + 1A 



This expression is valid only for sw/i) small so that 0o - sw/D. The term 
s/(w/2) represents the rate of decay due to holes recombining on the surface, 
s having the dimensions of velocity. For w » 1/too , the admittance becomes 
4w 2 qnp (iu/D) il2 , the same value as given in equation (4.12) for large o> and 
an area \w l . 

The conclusion from this appendix is that for x small, the effect of surface 
recombination is simply to modify the effective value of t and otherwise leave 
the theory of Section 4 unaltered. 

For very large values of x, it is necessary to consider higher terms in the 
sum and several values of r will be important. Under these conditions the 
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approximation is that, at x = 0, pi is independent of * and y may become a 
poor one, especially for forward currents, because the transverse currents to 
the edges will be important. Under these conditions the role of surface re- 
combination will give rise to patch effects of the sort discussed in Section 4. 

APPENDIX VI 

The Effect of Trapping Upon the Diffusion Process 

In this appendix we shall investigate the effect of the trapping of holes 
upon the impedance. We denote the density of mobile holes in the valence- 
bond band by p and the density of holes trapped in acceptors by p a . For 
thermal equilibrium at room temperature there will be an equilibrium ratio, 
called a, for p /p. For germanium a — 10 -4 and for silicon a — 0.1 to 0.2. 

We shall consider four processes which occur at rates (per particle per 
unit time) as follows: 

v r direct recombination of a hole with an electron (free or bound to a donor) 
vt trapping of a hole by an acceptor 
v ra recombination of a hole trapped on an acceptor 
v c excitation of a trapped hole into the valence-bond band. 

Under equilibrium conditions as many holes are being trapped (rate pvi) 
as are being excited (p a v e ): hence v t = av e . 

We shall study solutions of the customary form for the a-c. components: 

. iut—yx 

p\ = pwe 

. . fa ( JX 

P\a — Pla0 e 

These must satisfy the equations 

pi = Z>V"7>1 — (v t + V r )pi + Vepla 
Pla = Vlpl — (V e + V ra )pla 

These lead readily to the equation for y : 

Dy 2 = iiti + v r + v t — v e v t /(io) + v c + v ra ) = *'o> 

• L 1 + (,. + Zf + a, 2 ) J + Vr + " L 1 " (y. + O +"«7(*. + o J 

From this equation we can directly reach the important conclusion that 
the trapping process can never lead to a capacitative term larger than the 
resistive term. This result is obtained by analyzing the complex phase of 7, 
the admittance being proportional to y. In particular, we find that the real 
term in Dy 2 is always positive, as may be seen from inspection, so that the 
complex phase angle of 7 is less than 45°. 

The form reduces to a simple' expression if v e and v t are very large com- 



484 BELL SYSTEM TECHNICAL JOURNAL 

pared to v r , v ra and oj , a situation which insures local equilibrium between 
p and p a . Under these conditions we obtain 

Dy = Ki>[1 - a] + v r + ctv Ta 

Dividing by (1 + a) gives 

[d/{\ + a )W = [d p /(p + /> )] 7 2 = * + ^Vtr^ 

p ~r pa 

The interpretation is that the holes diffuse as if t ir diffusion constant were 
reduced by the fraction of the time p/(p + p a ) iey are free to move and 
recombine with a properly weighted average of v .nd v, a ■ 

APPENDIX VII 
Solutions of the Space Charge Equation 

We shall first show that the space charge equation (2.11) has a unique 
solution for the one dimensional case. For simplicity we write (2.11) in 
the form 

~ = sinh u - f(x) (A7.1) 

ax- 
to which it can be readily reduced. We shall deal with the case for which 
/ = / for.v<.r a (A7.2) 

/ = f b for x > Xb > Xa (A7.3) 

so that the interval (x a , x b ) is bounded by semi-infinite blocks of uniform 
semiconductor. We shall require that u be finite at x = ± <*> . This boundary 
condition requires that for large values of | * | 

11 = u a + A a e +yaX x -> - oo (A7.4) 

u = u b + A b e"' bX x -> + co (A7.5) 

where 

sinh u a = f a , sinh «& = f b 

7 „ = | (cosh u a ) 112 I , T 6 = I (cosh u b ) w I 

(If the opposite signs of the7's were present, the boundary conditions would 
not be satisfied.) The exponential solutions are valid for \u — Ua \ or 
| u _ Ub | « i. For larger values, however, solutions exist which are ob- 
tained by integrating (A7.1) to larger or smaller values of x. 

For these extended solutions the values of «(.v, -4 ) and u'(x, A a ) (= du/dx) 
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are monotonicaUy increasing functions of A a . This may be seen by con- 
sidering -v = x a . For A a sufficiently small, the value of u(x a , A a ) and u' (x a , A a ) 
are given simply by (A7.4). For larger values of A a , an exact integral will be 
required. It is evident, however, that 1 solutions of the form (A7.4) are 
related simply by translation for x < x a . Hence increasing A a is simply 
equivalent to integrating (A7.1) to larger values of x and it is evident that 
this increases u and u' mo otonically. It may be verified that for a sufficiently 
large A a the solution becc cs infinite at .v„ so that u(x a , A a ) u'(x a , A a ) both 
vary monotonicaUy and c tinuously from — °o to + » as A a varies from 
negative to positive valu We shall refer to this property of u(x a , A a ), 
u'(x a , A a ) as Pi . 

We next wish to show that u(xi , A a ), u'(xi , A a ) has the property Pi for 
values of X\ > x a . To prove this we note that if for any X\ , u(xi , A a ) and 
u'(xi , A a) are finite, the solution may be integrated somewhat further to 
obtain u(x-> , A a ), n'(x» , A a ) for x 2 > %\ . From equation (A7.1) it is evident 
that an increase in either u(xi , a) or u'(xi , a) will result in an increase in 
d 2 u/dx 2 in the interval Xi < x < x 2 so that u and u' at x% are monotonicaUy 
increasing functions of u and u' at X\ . Hence if u and u' at Xi have the 
property Pi , so do n and u' at .r 2 . By extending this argument we conclude 
that u and u' at any value of .v have the property Pi . (A rigorous proof 
can easily be completed along these lines provided that | /(*) | is finite.) 

Similarly it may be shown, starting from (A7.5), that u(x, Ab) is a mono- 
tonicaUy increasing function of At and u'(x, Ab) is a monotonicaUy decreasing 
function of Ab . 

In order to have a solution satisfying (A7.4) and (A7.5) we must have, 
for any selected point ;v, 

u(x, A a ) = u(x, A b ) (A7.6) 

u'(x, A a ) = u'(x, A b ) (A7.7) 

Now as the equation u(x, A a ) = u(x, Ab) varies from — oo to + =», u'(x, A a ) 
varies from — « to + °o and u'(x, Ab) varies from + co to — oo, monotoni- 
caUy and continuously. Hence there is one and only one solution of (A7.1) 
satisfying (A7.4) and (A7.5). 

In order to verify that the solutions discussed in Section 2 are correct for 
large and for small K, we show schematically in Fig. Al the solution for a 
representative K as a dashed line together with the curve u = i(o(y) = sinh - 
y. In terms of «o , equation (2.16) becomes 

— = — (sinh u — sinh « ). (A7.8) 

dy- K- 
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From the symmetry of the equation, it is evident that u must be an odd 
function of y and hence that the solution must pass through the origin. 
The boundary condition in this case will be that u -+ u a for y — ► ± » so that 
there will be no space charge far from the junction. We can conveniently use 
the origin as the point at which the solution from y = + co joins that from 
y = — oo ; from symmetry, this requires merely that n = when y = 0. 



U. AND U-o T 


u =si.nh'y ^ — -^Z^ 

/ S 

/ s 


/ / 

* / 

,' / 
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y 

df^ = ' [sLnhu-sinhuol 
dy 2 K 2L J 



Fig. Al— Behavior of the solution of Equation (2.16) or (A7.8). 



For large negative y, u = sinh -1 y and du/dy = 1/cosh m so that rfM/(/y 
is small. It is at once evident that, for large values of K, u must lie above « 
so that the integral 



c' J du 

(l/K 2 ) I (sinh U - sinh w ) dy = j- 



(A7.9) 



will be large enough to make the solution u(y) pass through the origin. If 
u — Uo> 2 over the region of largest difference, the space charge will be 
largely uncompensated and the solution will correspond to that used in 
equation (2.18). On the other hand, as K -> 0, the requirement that u(y) 
pass through the origin leads to the conclusion that u - Mo must be small for 
all values of y. The possibility that u oscillates about Mo need not be con- 
sidered since it may readily be seen that, if for any negative value of v, 
say yx , both u(yi) and «'()i) are less than « (yi) and u'(yi), then u(y) and 
u'iy) are progressively less than w (y) and u' (y) as y increases from y x to 0. 
Hence, if for negative y the u curve goes below the u curve, it cannot pass 
through the origin. 
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APPENDIX VIII 

List of Symbols 

(Numbers in parentheses refer to equations) 

a = (N d - N a )/x (2.14) 

.1 = admittance per unit area of junction (4.23) 

A p = component of A due to hole flow into w-region (4.12) (4.24) 

A n = component of A due to electron flow into />-region (4.25) 

At — component of A due to varying charge distribution in transition 

region 
A also used as a constant coefficient in various appendices 
b = ratio of electron mobility to hole mobility 
b = symbol for base in Sections 5 and 6 
B constant coefficient in various expansions in appendices 
c = symbol for collector in Section 6; a length in Appendix III 
C = capacity per unit area 
C„ , C p (4.25) (4.27) as for A n , A p 
C T (2.42) (2.45) (2.56) as for A T 

D = diffusion constant for holes (bD is the diffusion constant for electrons) 
e = 2.718... 
/ see Appendix 7 

g = rate of generation of hole-electron pairs per unit volume (3.1) 
G = conductance per unit area of junction 
G„ , G p as for ^4's 
*'= V-i" 
/ = current density 

In , I ,, = current densities due to electrons and holes (2.5) (2.6) (4.10) 
Ino.Ipa / P i (4.11) (4.12) (4.18) (4.19) 

/• , /n«j I P s saturation reverse current densities (4.11) (4.18) (4.21) 
I r see text with (4.35) 

/ = subscript in Section 3 for junction Fig. 5 equation (3.11) 
k = Holtzmarin's constant 
K = space charge parameter (2.17) 
L = length 
L a = Hi/a (2.15) 
L D = Debye length (2.12) 

L n , L p = diffusion lengths for electron in ^-region and holes in w-region (4.8) 
I. r = length required for potential increase of kT/q in region of constant 

field (4.32) Appendices II and IV 
L\ corresponds to a-c. held, Appendix IV 
n = density of electrons 



488 BELL SYSTEM TECHNICAL JOURNAL 

n n , n p = equilibrium densities of electrons in n- and />-regions 

p = density of holes 

Pn , p P = equilibrium densities of holes in //- and ^-regions 

p = d-c. component of non-equilibrium hole density (4.3) 

pi exp(ico/) = a-c. component of non-equilibrium hole density (4.3) 

p = total number per unit area of holes in specimen (2.35) 

q = electronic charge (q = \q\) 

Q = qP = total charge per unit area (2.39) 

r = recombination coefficient for holes and electrons (3.1) 

R = resistance of unit area 

R = resistance of unit area obtained by integrating conductivity (3.10), 

Appendix I 
Ri = effective series resistance, discussed in connection with (3.13) 
s = rate of recombination per unit area of surface per unit hole density, 

Appendix V 
S = susceptance per unit area (imaginary part of admittance) 
S pt 5«, S T as fork's. 
/ = time 

T = temperature in °K 
T = subscript for transition region 
u = $/kT (2.9), qty - tpi)/kT (2.32), Appendix VII 
vo and Vie iut = d-c. and a-c. components of voltage applied in forward direc- 
tion (4.2) 
W = width of space charge region in abrupt junction, Section 2.4 
w = half thickness of w-region or transistor base of Sections 5 and 6. 
■w = half width of square rod in Appendix V. 
x = coordinate perpendicular to plane of junction 
y, z = transverse coordinates, Appendix V 
y = reduced length (2.17), Appendix VII 
a = current gain factor in transistor (6.4) 
a = parameter in Appendix III and VI 
an = parameter in Appendix V 
/J t - = parameter in Appendix V 
7 = parameter in Appendices II, IV and VII 
e = symbol for emitter Section 6 
6,- = Pflu Appendix V 
k = dielectric constant 

H = mobility of a hole (6/i = mobility of electron) 
v = rates of recombination etc., Appendix VI 
p = charge density (2.1) 
a = conductivity 
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a i = conductivity of intrinsic material (4.15) 

<r„ = conductivity of H-region = qbnn n 

<t p = conductivity of /^-region = qnp p 

t = time 

t„ | t p = life times of electrons in /(-region and holes in w-region (3.2) (3.3) 

(4.7) 
Tr= relaxation time of transition region, Appendix IV 
tp, <p P , $n = Fermi level and quasi Fermi levels (2.2) (2.4) 
8<p = applied voltage across specimen in forward direction, Section 2.3, 

(4.2) 
X = sw/D in Appendix V 
\p = electrostatic potential (2.2) 
co = circular frequency of a-c. (4.2) 



